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The thermodynamical property of a small cluster including M Hubbard dimers, each of 
which is described by the two-site Hubbard model, has been discussed within the nonex- 
tensive statistics (NES). We have calculated temperature and magnetic-field dependences of 
the specific heat and susceptibility for M = 1,2, 3 and oo, assuming the relation between 
M and the entropic index q given by q = 1 + 1/M, which was previously derived by several 
methods. For relating the physical temperature T to the Lagrange multiplier j3, two methods 
have been adopted: T = l/k B /3 in the method A [Tsallis et al. Physica A 261 (1998) 534], 
and T = c q /k B /3 in the method B [Abe et al. Phys. Lett. A 281 (2001) 126], where k B 
denotes the Boltzman constant, c q — J^pf, an d Pi the probability distribution of the ith 
state. A comparison between the results calculated by the two methods suggests that the 
method B may be more appropriate than the method A for small-scale systems. 

§1. Introduction 

In the last several years, much study has been made with the use of nonextensive 
statistics (NES) which was initiated by Tsallis. 1 ' 4 ) Before discussing the NES, let's 
recall the basic feature of the Boltzman-Gibbs statistics (BGS) for a system with 
internal energy E and entropy S, which is immersed in a large reservoir with energy 
Eq and entropy Sq. The temperature of the small system T is the same as that of 
the reservoir To where T = 5E/5S and To = 5Eq/5Sq. If we consider the number 
of possible microscopic states of Q(Eq) in the reservoir, its entropy is given by 
So = ks \nQ(Eo) where ks denotes the Boltzman constant. The probability of 
finding the small system with the energy E is given by p(E) = Q{Eq — E)/f2(E) ~ 
exp(— E/ksT) with E <C Eq. When the physical quantity Q of a system containing 
N particles is expressed by Q oc iV 7 , it is classified into two groups in the BGS: 
intensive (7 = 0) or extensive one (7 = 1). The temperature and energy are typical 
intensive and extensive quantities, respectively. This is not the case in the NES, as 
will be shown below. 

In the NES, on the contrary, the temperature of a nanosystem which is in contact 
with the reservoir, is expected to fluctuate around the temperature of the reservoir 
To because of the smallness of nanosystems and their quasi-thermodynamical equi- 
librium states with the reservoir. Then the BGS distribution mentioned above has to 
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be averaged over the fluctuating temperature. This idea has been expressed by 5 ) 6 ) 7 ) 

poo 

p(E) = / d/3e-f> E f B (0) 
Jo 

= [1 - (1 - q)0 o E]^ = exp q (-0 o E), (1-1) 



with 



9 = l + £. (1-2) 

N 

1 / N \~* -N-i ( N/3 



0o = -4- = f°° d0f(0) = E(0), (1-4) 
kbTq Jo 

2 E{0 2 ) - E{0) 2 



N E(0) 2 
where exp 9 (x) denotes the ^-exponential function defined by 



(1-5) 



exp g (x) = [1 + (1 - q)x] !-« , for 1 + (1 - q)x > 

= 0, otherwise (1-6) 

q expresses the entropic index, E(Q) the expectation value of Q averaged over the r 
(or x 2 ) distribution function f B (0), 0o the average of the fluctuating and 2/N its 
dispersion. The r distribution is emerging from the sum of squares of N Gaussian 
random variables, related discussion being given in Sec. 3. 

The functional form of the probability distribution p(E) expressed by Eq. (1.1) 
was originally derived by the maximum entropy method 1 ) 2 ) with the generalized 
entropy given by 

S q = k B C^fl' 1 ) = -k B £ Pi Kfa), (1-7) 

where pi [= p(ei)] denotes the probability distribution for the energy in the system, 
and ln„(x) [= {x l ~ q — 1)/(1 — q)\ the g-logarithmic function, the inverse of the q- 
cxponential function. The important consequence of the NES is that entropy and 
energy are not proportional to N in nanosystems. 

In our previous papers 8 )' 9 ) (referred to as I and II, respectively), we have applied 
the NES to the Hubbard model, which is one of the most important models in solid- 
state physics. Thermodynamical properties of canonical 8 )' 9 ) and grand-canonical 
ensembles 9 ) of a Hubbard dimer described by the two-site Hubbard model have been 
calculated within the NES. It has been shown that the specific heat and susceptibility 
calculated by the NES may be significantly different from those calculated by the 
BGS when the entropic index q departs from unity, the NES with q = 1 reducing to 
the BGS. 
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We will consider in the present paper, a nanocluster containing multiple (M) 
Hubbard dimers in order to discuss the M dependence of their thermodynamical 
properties. We have assumed the M — q relation: 

* = 1 + W (1 " 8) 

which is derived from Eq. (1.2) with N = 2M. 

The paper is organized as follows. After discussing the adopted model and 
calculation method, we will present in §2, numerical calculations of temperature and 
magnetic-field dependences of thermodynamical quantities for various M values. The 
final §3 is devoted to discussion and conclusion. 

§2. Nonextensive statistics for Hubbard dimers 

2.1. Adopted model and calculation method 

We have adopted canonical ensembles of a small cluster containing M Hub- 
bard dimers, each of which is described by the two-site Hubbard model. Interdimer 
interactions are assumed to be negligibly small. The Hamiltonian is given by 

M 

H = Y,Hf\ (2-1) 

2 2 
a j=l j=l 

(1,2 €£) (2-2) 

where denotes the Hamiltonian for the £th dimer, rij a = Oj a aj a , a,j a an annihi- 
lation operator of an electron with spin a on a site j (e £), t the hopping integral, 
U the intraatomic interaction, \xb the Bohr magneton and B an applied magnetic 
field. Six eigenvalues of are given by 

e t = 0, 2p B B, -2p B B, U, | + A, | - A, for % = 1 - 6 (2-3) 

where A = \/U 2 /4: + At 2 . 10 ^ 11 ^ The number of eigenstates of the total Hamiltonian 
H is 6 M . 

The entropy S q in the Tsallis NES is defined by 1 -* 2 ) 

Tr(p«)-1 



Here p q stands for the generalized canonical density matrix, whose explicit form will 
be determined shortly [Eq. (2.7)]. We impose the two constraints given by 

Tr (p q ) = 1, (2-5) 
-T^W = <H>q - ^ (2 ' 6) 
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where the normalized formalism is adopted. 2 ) The variational condition for the 
entropy with the two constraints given by Eqs. (2.5) and (2.6) yields 



1 

Pq = -=-exp, 

i 

with 

X q = Tr ( exp 



P 



( H - E, 



q) 



-)(H-E q ) 



(2-7) 



(2-8) 



c q = Tr{pl)=X 1 -\ (2-9) 

where exp g [x] expresses the (/-exponential function defined by Eq. (1.6) and (3 is a 
Lagrange multiplier: 

< 2 ' 10 » 

Specific heat and suscptibility have been calculated in I and II. 

For relating the physical temperature T to the Lagrange multiplier j3, we have 
adopted the two methods A and B, given by 9 ) 

T=-^—, (method A) (2-11) 

KBP 

= p-. (method B) (2-12) 

kbp 

The method A proposed in Ref. 2 is the same as the extensive BGS. The method B 
is introduced so as to satisfy the zeroth. law of thermodynamical principles and the 
generalized Legendre transformations. 12 ) It has been demonstrated that the negative 
specific heat of a classical gas model which is realized in the method A, 13 ) is remedied 
in the method B. 12 ) A difference between the two methods does not matter as far as 
we consider only the non-exponential distribution in the NES. It yields, however, a 
significant difference in the temperature dependence of thermodynamical quantities, 
as will be shown in the following subsection. 

2.2. Numerical calculations 
2.2.1. Temperature dependence 

In order to study how thermodynamical quantities of a cluster with Hubbard 
dimers depend on its size M, we have made numerical calculations, assuming the 
M — q relation given by Eq. (1.8), where results for M = oo correspond to those 
of the BGS ((7 = 1). Figures l(a)-l(d) show the results for U/t = 5. The specific 
heat and susceptibility shown in Figs. 1(a) and 1(b), have been calculated by the 
method A with q = 2.0, 1.5, and 1.333 for M = 1, 2 and 3, respectively. Figures 
1(c) and 1(d) express C q and Xqi respectively, calculated by the method B. We note 
that physical quantities in a small cluster with M ~ 1 — 3 are rather different from 
those of bulk-like systems with M = oo, although properties of clusters gradually 
approach those of bulk with increasing M. 

We note in Figs. l(a)-l(d) that by varying M, the maximum values of the 
specific heat (C*) and the susceptibility (x* q ) and corresponding temperatures of Tq 
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and T* are changed. Figure 2(a) shows Tq and T*, and Fig. 2(b) depicts C* and 
Xq, which are plotted against 1/M: solid and dashed lines denote results calculated 
by the methods A and B, respectively. It is shown in Fig. 2(a) that with increasing 
1/M, T* calculated by the method A is much increased than that calculated by the 
method B. We note also that with increasing 1/M, T c of the method B is increased 
while that of the method A is decreased. Figure 2(b) shows that C* in the method 
A is smaller than that in the method B, whereas x* q 1X1 t ne method A is the same as 
that in the method B. 

2.2.2. Magnetic-field dependence 

From the B dependence of the six eigenvalues of [Eq. (2.3)], we note the 
crossing of the eigenvalues of £3 and e§ at the critical filed: 

[u* U 

VBB c =y{— + * 2 - T , (2-13) 

leading to fiBB c /t = 0.351 for U/t = 5. For B < B c (B > B c ), e 6 (e 3 ) is the ground 
state. At B ~ B c the magnetization m q is rapidly increased as shown in Figs. 3(a) 
and 3(b) for ksT/t = 1.0 and 0.1, respectively: the transition at lower temperature 
is more evident than at higher temperature. This level crossing also yields a peak in 
Xq [Figs. 3(c) and 3(d)] and a dip in C q [Figs. 3(e) and 3(f)]. It is interesting that the 
peak of Xq f° r Q = 1-5 i n the NES is more significant than that in the BGS whereas 
that of C q of the former is broader than that of the latter. When the temperature 
becomes higher, these peak structures become less evident. Similar phenomenon 
in the field-dependent specific heat and susceptibility have been pointed out in the 
Heisenberg model within the BGS. 14 ) 

Figure 3(a) and 3(b) remind us the quantum tunneling of magnetization observed 
in magnetic molecular clusters such as Mn4, Mnl2 and Fe8. 15 ) It originates from 
the level crossing of magnetic molecules which are parallel and anti-parallel to the 
easy axis when a magnetic field is applied. 

§3. Discussions and conclusions 

The N — q relation given by q = 1 + 2/N [Eq. (1.2)] has been derived from 
the average of the BGS partition function over the T distribution f B given by Eq. 
(1.3). By using the large-deviation approximation, Touchette 16 ) has obtained the 
alternative distribution function / T (/3) given by 

Solid and dashed curves in Fig. 4 express the f B - and / T -distribution functions, 
respectively, for various iV values. For N — > 00, both reduce to the delta-function 
densities, and for a large iV = 100, both distribution functions lead to similar results. 
For a small ./V (< 10), however, there is a clear difference between the two distribution 
functions. We note that a change of variable (3 — > in f T yields the distribution 
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function similar to f B . It should be noted that f T cannot lead to the ^-exponential 
function which plays a crucial role in the NES. For a large e, f T leads to the stretched 
exponential form of p(e) ~ e c ^ while f B yields the power form of p(e) ~ e~^ =1 . This 
issue of f B versus f T is related to the super statistics, which is currently studied with 
much interest. 17 ) 

Numerical calculations presented in the preceding section have shown that al- 
though results calculated by the two methods A and B are qualitatively similar, there 
are some quantitative difference, as previously obtained in I and II. When we calcu- 
late the Curie constant r q of the susceptibility defined by X g (T) = (fi B / kB)[r q {T) /T], 
the method A leads to anomalously large Curie constant compared to that of the 
method B. 9 ) This agrees with the results for free spins 9 ) 18 ) and for spin dimers de- 
scribed by the Heisenberg model. 22 ) A comparison between Eqs. (1.1) and (2.7) 
yield the average temperature < T > given by 

, * ^ Po = - (3-2) 

k B <T > c q 

which is consistent with the method B. These results suggest that the method B 
is more appropriate than the method A. This is consistent with recent theoretical 
analyses 19 ) 20 ) [for a relevant discussion, see also Ref. 21]. 

In summary, within the framework of the NES, thermodynamical properties have 
been discussed of a cluster including M dimers, each of which is described by the two- 
site Hubbard model. We have demonstrated that the thermodynamical properties 
of small-scale systems are rather different from those of bulk systems. Owing to 
recent progress in atomic engineering, it is possible to synthesize molecules containing 
relatively small numbers of magnetic atoms with the use of various methods (for 
reviews, see Refs. 21-23). Theoretical and experimental studies on nanoclusters 
with changing M could clarify a link between the behavior of the low-dimensional 
infinite systems and nanoscale finite-size systems. The unsettled issues on T — (3 
and the N — q relations in the current NES are expected to be resolved by future 
theoretical and experimental studies on nanosystems, which are expected to be one 
of ideal systems for a study on the NES. 
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Fig. 1. The temperature dependences of (a) specific heat C q and (b) susceptibility \q (P er dimer) 
of Hubbard dimers for U/t = 5 calculated by the method A, and those of (c) specific heat C q 
and (d) susceptibility \q calculated by the method B, with M = 1 (bold solid curves), M — 2 
(chain curves) , M = 3 (dashed curves) and M = oo (solid curves) . 



Fig. 2. (a) 1/M dependence of the temperatures of Tq (circles) and T* (squares) where C q and 
Xq have the maximum values, respectively, (b) 1/M dependence of the maximum values of C* q 
(circles) and Xq (squares). Solid and dashed lines denote the results calculated by the methods 
A and B, respectively: T* calculated by the method A shown in (a) is divided by a factor of 
five. 



Fig. 3. The magnetic-filed dependence of (a) the magnetization m q for fcsT/t = 1.0 and (b) 
ksT/t = 0.1, (c) the susceptibility for ksT/t = 1.0 and (d) ksT/t = 0.1, (e) the specific heat 
X q for k B T/t = 1.0 and (f) k B T/t = 0.1, with U/t = 5 for M = 2 calculated by the method A 
(solid curves) and B (dashed curves) and for M = oo (chain curves). 



Fig. 4. The distributions of f B (f3) (solid curves) and f T (f3) (dashed curves) as a function of f3 (see 
text). 



